Transition from isotropic to anisotropic beam profiles in a linear focusing channel. 
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This paper examines the transition from isotropic to anisotropic beam profiles in a linear focusing 
channel. Considering a high-intensity ion beam in space-charge dominated regime and large mis- 
matched RMS beam size initially, observe a fast anisotropy situation of the beam, characterized for 
a transition of the transversal section round to elliptical with a coupling of transversal emittance 
driven for instabilities of nonlinear space-charge forces. Space-charge interactions in high-intensity 
linear accelerator can lead to equipartitioning of energy between the degrees of freedom. The 
anisotropization phenomena suggest a kind of route to equipartition. In order to understand the 
initial dynamical behavior of an anisotropic beam, in particular, to study possible mechanisms of 
equipartition connected with phase space we have to know how we can compute the variables (vol- 
ume, area of surface, and area projected) that characterize the anisotropic beam in phase space. 
The beam began in a nonequilibrium state evolved toward a metaequilibrium in which the particle 
orbits filled an invariant measure of phase space. We shall call this subspace of the phase space 
" ^ ", where ^ is the ratio of oscillations energies in the x and y directions. It is well-know [s^ that 
an isolated difference ressonance of the form luj^ — mujy = 0, where uox is single particle frequency 
along of a;-direction and ujy is single particle frequency along of y-direction, lead to the motion 
bounded in both directions and " ^ "remains unchanged. The purpose of this paper is to propose 
one definiton of the anisotropic equipartition [sj . Anisotropic equipartition corresponds to a phase 
space density uniform on the surface invariant of the f, a version of the ergodic hypothesis where 
the ^ invariant play the role of the conserved energy In the state of anisotropic equipartition, 
the beam temperature is stationary, the entropy grows in the cascade form, there is a coupling of 
transversal emittance, the beam develops an elliptical shape with a increase in its size along one 
direction and there is halo formation along one direction preferential. 

PACS numbers: 41.75.-i, 41.85.Ja, 29.27.Bd, 52.59.Sa, 52.25.Kn, 45.50.Dd 



I. INTRODUCTION 

In space-charge dominated beams the nonlinear space- 
charge forces produce a filamentation pattern, which re- 
sults in a 2-component beam consisting of an inner core 
and an outer halo P, S S I3- When this core is mis- 
matched in a uniform linear focusing channel, the en- 
velope oscillates and the particles, represented by single 
test particles, oscillate about and through the core. This 
mechanism is called particle-core model, considering a 
test particles initially located outside the core. The test 
particles execute betatron oscillations under the influence 
of space-charge field induced by the oscillating core, and 
exhibit various nonlinear behaviors including paramet- 
ric resonance. In a number of numerical analyses and 
macro-particle simulations, it has been revealed that the 
2:1 parametric resonance between test particle oscillation 
and breathingcore oscillation is a main cause of halo for- 
mation @, i, 0, i, i, 0, [H • 

Space-charge induced coupling between different de- 
grees of freedom can be responsible for emittance growth 
or transfer of emittance from one phase plane to another. 
The underlying instability mechanism is coherent if it 
depends primarily on the electric field due to collective 
motion [l^, [13] , in contrast with an incoherent that can 
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be described as single-particle effect. In an analytical 
single-particle analyses Montague [12 pointed out that 
the space-charge driven fourth order difference resonance 
may lead to emittance coupling. The theory of space- 
charge coupling resonances has been studied most thor- 
oughly by Hofmann using Vlasov equation [l5l | . Intrinsic 
coupling resonances are an important consideration in 
linacs and may also be of interest in rings when the tune 
separation is small [16]. These resonances are accompa- 
nied by halo formation in the receiving energy. In the 
presence of internal energy anisotropic between different 
degrees of freedom, initially small space-charge coupling 
terms can grow exponentially due to collective instability. 
The coupling resonaces driven by the beam space-charge 
fields depending only on the relative emittances and av- 
erage focusing strengths. In anisotropic beams, the emit- 
tance and/or external focusing force strength are differ- 
ent in the two transverse directions. Ikegami [2lj deal 
with effects of anisotropic of beam cores on halo dynam- 
ics. 

Most of the halo studies so far have considered round 
beams with axisymmetric focusing. Some news aspects 
caused by anisotropy demonstrate an influence of the 
mismatch on halo size 13, 26] . The mismatch oscillations 
can drive particles into a halo as a result of resonant inter- 
action of these particles with the mismatch mode [6^ . So 
far only second order round beams have been considered 
as possible mismatch modes: the influence of anisotropic 
on second and higher order mismatch modes is expected 



2 



to be an important factor in halo formation. 

Previous work has been done in a nonhnear analyses of 
the beam transport considering nonaxisymmetric pertur- 
bations [23 |. In which it is shown that large-amplitude 
breathing oscillating of an initially round beam couple 
nonlinearly to quadrupole-like oscillations, such that the 
excess energy initially constrained to the axisymmetric 
breathing oscillating is allowed to flow back and forth 
between breathing and quadrupole-like oscillations. In 
this case, the beam develops an elliptical shape with a in- 
crease in its size along one direction as the beam is trans- 
ported. This is a higly nonlinear phenomenon that occurs 
for large mismatch amplitudes on the order of 100% [ioj . 

This papers [23] deals with the coupled motion between 
the two tranverse coordinates of a particles beam arising 
from the space-charge forces, and, in particular, with the 
effects of the beam coupling which have ratio r^jry — 1, 

and Ty are the envelope semi-axes-rms. A beam with 
nonuniform charge distribution always gives rise to cou- 
pled motion. However, it is only when the ratio r^/ry = 1 
and the large beam size-rms mismatched on the order of 
100% (4flJ that the coupling can produce an observable 
effect in the beam as a whole. This effect arises from 
both, a beating in amplitude between the two coordinate 
directions for the single-particle motion, and from the 
coupling between oscillations modes beam, resulting in 
growth and transfer of emittance from one phase plane 
to another, in the beam develops an elliptical shape and 
therefore in the transition from isotropic to anisotropic 
beam profiles. 

Nonlinear space-charge forces can also lead to equipar- 
titioning of energy between degrees of freedom. The ques- 
tion is how a system of collisionless particles cou pled by 
long-range space-charge forces will equipartition. [3l| . In 
the presence of nonlinear coupling mechanisms, the beam 
can be predicted to equipatition, or reach a state where 
velocity spreads in two directions are equal. In space- 
charge dominated beams. Coulomb collisions are infre- 
quent to account for the energy transfer, whereas space 
charge waves have been shown to be a possible coupling 
canditate [1^, [l^, H^. Equipartitioning of anisotropic 
beam involves nonlinear energy transfer and evolution 
towards a metaequilibrium state, as a consequence of 
resonant phase mixing [i^, [gO, EH. Strictly speaking, 
resonant phase mixing is a reversible process in which it 
is governed by Vlasov's equation. However, an essential 
question for the accelerator designer is wheter this pro- 
cess in operationally reversible. While it may be possible 
in principle to compensate operationally against phase 
space dilution, this compensation must be completed be- 
fore any mixing has smeared a significant number of par- 
ticles through global regions of phase space. It arises re- 
garding any process for manipulating a beam with space- 
charge [25|, be this changing the beam's tranverse ge- 
ometry (round-to-elliptical or elliptical-to-round trans- 
formations [l^iHJ)- We are working toward quantifying 
the relationship between the anisotropy of the beam and 
the equipartition. The equipartition of beam is driven 



for anisotropics processes. It is, therefore, to develop an 
improved understanding of fundamental collective stabil- 
ity properties, including the case where a large temper- 
ature anisotropy can drive electrostatic Harrys-type [s^ 
and/or eletromagnetic Weibel-type [sj instabilities, fa- 
miliar in the study of electrically neutral plasmas [55|. 
In plasmas with strongly anisotropic distribution func- 
tions, collective instabilities may develop if there is sufh- 
cient coupling between the degrees of fredom. Previous 
studies have mostly focused on the electrostatic Harris- 
type anisotropic-driven instability for beams [56| . It has 
been shown that a fast, electrostatic instability devel- 
ops, and satures nonlinearly, for sufficiently temperature 
anisotropic. 

The term equipartition broadly refers to the er- 
godic property of multi-dimensional Hamiltonian sys- 
tems, which tend to distribute uniformly over the phase 
space surface of constant energy. The conservation of 
energy plays the fundamental role in classical equilib- 
rium thermodynamics. The term "turbulent equiparti- 
ton "was introduced by Yankov [ss'l in order to describe 
the turbulent relaxed state, in which the system assumes 
a uniform distribution on the surface of constant invari- 
ants respected by turbulence. In plasmas physics, the 
best known example of a turbulent equipartition is the 
quasilinear plateau of the distribution function caused 
by the nonlinear Landau damping of plasmas waves. In 
a toroidal turbulent plasma, the relevant invariants are 
given by the frozen-in law (in the fluid limit) or the adia- 
batic invariant and the Liouville theorem (in the collision- 
less limit), both limits derivable from the more general 
Poincare invariant. The plasma mixing by low- frequency 
electrostatic modes in a Tokamak, subject to these con- 
servative laws, results in the inhomogeneous density and 
temperature profiles peaked at the center even in the ab- 
sence of particles and energy fluxes, thus presenting the 
underlying mechanism of the pinch effect and the pro- 
files consistency in Tokamak [s^, HI]- The purpose of 
this paper is to propose one definiton of the anisotropic 
equipartition. Anisotropic equipartition corresponds to 
a phase space density uniform on the surface invariant 
of the (^ = [ryex)"^ /{rx^y)"^): where ^ is the ratio of en- 
ergy oscillations in the x and y directions, a version of 
the ergodic hypothesis where the ^ invariant play the role 
of the conserved energy [6j | . In the state of anisotropic 
equipartition, the temperature is stationary, the entropy 
grows in the cascade form, there is a coupling of transver- 
sal emittance, the beam develops an elliptical shape with 
a increasing in its size along one direction and there is 
halo formation along one direction preferential. 

This paper is organized as follow. In Sec. [Til the models 
equations are derived. In Sec. Ill Al we examine the tran- 
sition from isotropic to anisotropic beam profiles in a lin- 
ear focusing channel and the transition of the transver- 
sal section from round to elliptical, with a coupling of 
transversal beam emittance. In Sec. lIIII we study possible 
mechanisms of equipartition connected with phase space 
and compute the variables (volume, area of surface, and 
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area projected) that characterize the anisotropic beam in 
phase space. FinaUy in Sec. IIVI we discuss impHcations 
and possibles extensions of our results. 

II. THE MODEL EQUATIONS 

We consider an axially long unbunched beam of ions 
of charge q and mass m propagating with average ax- 
ial velocity (3bcez (c is the speed of light in vacuo and 
relativistic factor 7;, — — /3^) along an uniform lin- 

ear focusing channel, self-field interactions are electro- 
static. The beam is assumed to have an elliptical cross 
section centered at a; = = y and vanishing canoni- 
cal angular momentum Pg = {xy — yx ) = 0, where x 
and y are the positions of the beam particles. We con- 
sider nonuiform density beam in space-charge dominated 
regime and Cx = emittance, initially. The general 
property of space-charge dominated beam behaviour is 
that a beam with an initial nonlinear profile tends to be 
more uniform and this process is associated with strong 
emittance growth and the appearance of beam halo. 



As demonstrated by Sacherer [11] and Lapostolle [i^ 
enevelope equations for a continuous beam are not re- 
stricted to uniformly charged beams, but are equally 
valid for any charge distribution with elliptical sym- 
metry, provided the beam boundary and emittance 
are defined by rms (root-mean-square) values. Thus, 
we consider the parabolic density beam (n;, = 
2Affc/7rr^ry [l — x^/r^ — y^/r^] ) where iVf, is the axial 
line density, Vx = \jQ{x'^) and Vy — y^6{y^) are ellipsis 
semi-axes rms. A main point is that for parabolic den- 
sity distribution the fourth order space-charge potential 
driving the coupling is already present in the initial dis- 
tribution, hence emittance exchange and beam develops 
an elliptical shape immediately. 

For a parabolic density n{x,y), Poisson's equation 
V'^cj) — —qn/eo, So is the permittivity of free space, pro- 
vides the basis for obtaining the space-charge field com- 
ponent (assuming paraxial approximation). The density 
is assumed to be zero outside of the ellipse. The solution 
has been given by Lapostolle The electrostatic po- 
tential (j) is given by : 
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A = (a;2 - y2 - A2)2 -f Ax'^y^ and A='= = 



VA±(a;2-y2zpA2). 

The transverse orbit x{s) of a beam particle satisfy the 
paraxial equation of motion 
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with an analogous equation for orbit y{s). Here, s is the 
axial coordinate of a beam, primes denote derivates with 
respect to s and kq is represented constant focusing force. 

The envelope of the beam is an elliptical cross-section 
with rms radii rj (henceforth, j ranges over both x and 



y) that obey the rms-ivTV^ envelope equations 



rj + Kgrj 



2K 
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Here, K = q Nij/tt e()j'^f3^mc is the dimensionless per- 
veance of the beam, is rms-emittance of the beam 
along the j-plane. 

The Ej = V — {jj')'^ can analytical been cal- 
culated following a model proposed to Lapostolle et 
al. [A2\ for nonlinear space-charge forces, and a nonlinear 
analysis proposed to Pakter et al. (2^ for large mismatch 
beam. According to Lapostolle et al. |42l | nonlinear space- 
charge forces cause a change in the momentum compo- 
nents, which is equal to the product of the force and the 
time over which the force acts. The force depends on the 
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spatial distribution of the particles, and particle coordi- 
nates on which the force acts. In general, these changes 
in the momentum components change the phase-space 
distribution of the particles. Thus the particles experi- 
ence a space-charge impulse, but do not propagate far 
enough for their positions to change appreciably. For a 
beam with large mismatch amplitudes, Pakter et al. (2^ 
has demonstrated that mismatch is like an impulse so 
fast that particle's positions not change, only momentum 
suffers a discontinuous variation. Mismatched oscillato- 
ries modes beam oscillate between a maximum value and 
a minimum value around matched envelope with a given 
periodicity. To label different mismatched oscillations, its 
define a mismatch amplitude as A = ^max/ipo > 1 where 
Tpmax is maximum oscillation amplitude of breathing and 
quadrupole modes, and i/iq is matched oscillation ampli- 
tude. Mismatched oscillations are excess energy given to 
the beam. In particular, there is a threshold mismatched 
amplitude At^ above which effective energy exchange be- 
tween breathing and quadrupole modes takes place [2^ . 

From electrostatic potential and mismatch amplitudes 
on the order of 100%, the tranvcrse momentum impulse 
can be calculated. This results in a new phase-space 
distribution and new rms emittance j43l |44 |. For ex- 
ample, in the x plane; the change in the momentum 
component is ISp^ = qExA/vjy, where = d(t)™/dx 
is electric field, Wfc is the beam velocity. The impulse 
can also be expressed as change in the divergence angle, 
given non-relativistically in the paraxial aproximation by 
Aa; = qExA/rribv'^. If the second moments of the parti- 
cle distribution can be evaluted from the expression for 
X and X , the rms emittance can be obtained. Know- 



ing that divergence is a; = x/r^ + qE^Ai 
emittance for parabolic density yields [231 



the rms 



1/2 



(1 + ^; 
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the result is easily transformed to the y plane interchang- 
ing rx and ry. It observes that emittance depends of 
perveance beam K, mismatched amplitude of oscillato- 
ries modes beam A and ratio of the beam semi-axes. The 
first term corresponds to the filamention effect caused by 
the fourth order term in the electrostatic potential. The 
second term comes from the coupling, i.e the dependence 
of the X component of the potential on the y coordinate, 
which produces the spreading of the initial filamention. 
The last term is a cross term between the filamention 
term and the coupling term. 

KV distribution has frequently been taken as a theo- 
retical basis of particle-core model because space-charge 
forces are linear. However, there is no doubt that real- 
istic intense beams contain a fully nonlinear nature. It 
may thus be reasonable to try constructing an alternate 
particle-core algorithm with a nonlinear core potential 
even if the model is only aproximation. From this point 



of view, we introduce here a parabolic core under a sim- 
plifying assumption . In this particle-core model, the core 
is described by rms envelope equation ([U, and the halo 
particles are modeled using test particles that subject 
to the external force and the time-dependent nonlinear 
space-charge force associated with the parabolic core. We 
assume that the parabolic-type density is roughly main- 
tained even for a mismatched beam. Thus the spatial dis- 
tribution is assumed to remain unchanged as the beam 
propagate. The test particles do not affect the motion 
of the core [i^ and they are described by equation ([3|) 
taking for inside beam, and (j)out for outside beam. 
Thus the theory is not self-consistent. 



A. Transition from isotropic to anisotropic beam 

It is easy to verify that there is a particular solution 
of the envelope equations (|4]) for which rj{s) = t'^q = 
[{K + {K^ + 4^2^2)1/2-) /2k2] 1/2^ ^ijgre ^ ^ ^^/^^ 
tio emittance. This corresponds to the so called matched 
solution for which a circular beam of radius r^Q pre- 
serves its shape throughout the transport along the fo- 
cusing channel. Then, we transform the equations to 
a dimensionless form introducing the following dimen- 
sionless variables an parameters: r = kqs for the inde- 
pendent variable, = yW/e^^'x and fy — ^ K^/tyry 
for envelope beam, x = ^ K^/eyX and y = ^Ko/cy?/ 
for test-particle, and K = K/eyKo for the scaled space- 
charge perveance. In addition, we introduce the follow- 
ing anisotropy variables: the ratio emittance 77, ratio of 
the envelope beam x = fx/^y and the mismatch factor 

= rx/no = Ty/rto. 

We launch the beam with K = 3, kq = 1, i' = 2.4, A = 
2.0, rjQ = vrto, Txo = ryo and rj = 1 initially, and inte- 
grate firstly the envelope equations (|4]) up to s = 50.0. 
The corresponding evolution of the rms-emittance ([5]) 
is used for ej in Eq.(|3]). It is convenient to introduce 
new canonical variables defined as Xg = (r^ + f y)/2 and 
Xa — (tx — Ty)/2 to analyze oscillations modes beam. 
Note that Xg are oscillations where rx{s) and ry{s) os- 
cillate in phase : breathing modes. And Xa are oscilla- 
tions where rx{s) and ry{s) oscillate with opposite phase 
: quadrupole modes [51|. 

In Fig. [T] the beam develops an elliptical shape with a 
increase in its size along of x-direction [2^. This effect 
is caused for large beam size-rms mismatched its that 
coupled the oscillations modes beam, perturbing nonlin- 
ear space-charge force [48| . Therefore perturbation in- 
duction coupling between different degrees of freedom 
and core-core resonances [l^, [s^. Real core-core reso- 
nances jisj are observed for rms mismatched beam ful- 
filling internal resonance conditions between the planes. 
Beam size-rms mismatched is excess energy given to the 
system. In general, this excess energy appear as oscil- 
lation energies in the degrees of freedom of the system. 
Space-charge couple some degrees of freedom causing res- 
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FIG. 1; Evolution of the envelopes obtained by direct inte- 
gration of the envelope equations Q (top). Vx and ry are 
represented by red and blue lines, respectively. Evolution of 
the oscillations modes beam (bottom) . Xs and Xa are repre- 
sented by red and blue lines, respectively. 



onance between them. Coupling resonance leading to an 
exchange of energy between the relevant degrees of free- 
dom. In Fig. [1] (bottom) there is an increase of oscilla- 
tion amplitude of modes featuring a resonance between 
breathing and quadrupole modes. To analyse this reso- 
nance we computed dimensionless frequencies of breath- 
ing and quadrupole modes using Fourier analyses these 
modes. The frequency associated with the breathing or 
quadrupole mode corresponds to the maximum in the 
Fourier transform. We used FFTW (The Fastest Fourier 
Transform in the West) to compute numerically the fre- 
quencies. Features this method are found on the web 
page [57] and papers 58]. Breathing mode frequency 
is LUx, = 2.01334 and quadrupole mode frequency is 
ujx^ — 2.01334. Note that both modes has the same 
frequency. Therefore breathing mode and quadrupole 
mode are resonance. This core-core resonance together 
with single-particle resonances are making the beam to 
develop an elliptical shape with a increase in its size 
along of x-direction [l^ . It should be noted that in the 
anisotropic case (x 7^ 1 and 77 7^ 1) both the breathing 



and quadrupole modes have quadrupolar symmetry. 

The emittance evolution given by equation ^ is shown 
in Fig. [2] The initial transversal emittance is equal 
^xo — Eyo — 0.22360 [23]. It is observed emittance cou- 
pling caused for space charge driven nonlinear difference 
single particle resonance [12, HI]- This coupling is 
characterized by the emittance exchange between the di- 
rections. Space-charge alters the net force seen by the 
individual particles in a way that is nonlinear and de- 
pendent on the density distribution of the beam. Emit- 
tance exchange requires resonant coupling, which can 
take place only if an intrinsic resonace relationship is 
fulfilled. A simplified approach would be a difference 
single particle resonance condition lujx — mujy = which 
was suggested in Ref. [3l|. ujx is single particle dimen- 
sionless frequencies along of x-direction and ujy is single 
particle dimensionless frequencies along of y-direction, 
I and m are integer numbers. The motion is always 
bounded in both directions . This emittance coupling is 
a manifestation of the exchange of energy from one to the 
other directions which is familiar in the linear coupling 
l = m=l. We used FFTW (The Fastest Fourier Trans- 
form in the West) to compute numerically the frequencies 
ujx and ojy of 2500 test particles. We launch the test par- 
ticles along the x- and y-axes of beam in specified region 
inside the beam (between O.Olrx and O.lvx spaced by 
0.0004ra; along x-direction and between O.Olrj, and O.Try 
spaced by 0.0004rj, along y-direction) . These 2500 test 
particles has uJx = 1.05461 and 2455 test particles has 
LUy = 1.05461. Thus 2455 test particles obey resonance 
condition lu^ — iUy = 0. The results are presented by the 
histogram in Fig. 3) This number large of test particles 
in resonance cause emittance coupling. The emittance 
growth is induced by nonlinear space-charge forces, 
increases and Cj, decreases with increasing r^/ry but ey 
increases and ex decreases with decreasing rx/ry as ob- 
served in Fig. [5] (bottom). The emittance growth is larger 
in the plane that has the larger semi-axes length. That 
the emittance growth increases as the semiaxes length 
increases may seem surprising for an effect that arises 
from space-charge force, which increases as the beam 
size becomes smaller, rather than larger and vice versa. 
i^x + ei/)/2 shows a much smoother but similarly pro- 
nounced anisotropic variable x response. 

An important effect of space charge is the tendency to 
induce waves the beam, a collective effect. These waves 
are characterized by the plasma frequency, which in turn 
relates to the perveance K defined as the ratio of the 
average space charge force to the external focusing force 
at the beam edge. The effect these waves is observed in 
the asymtotic evolution of the envelope and emittance 
shown in Fig. [31 The exponential growth of the envelope 
and the emittance exchange is characteristic by insta- 
bility tilting mode of a space-charge-dominated regime 
beam (l^ [la ] . Space-charge in connection with nonlin- 
ear dynamics lead to this phenomenon. The parabolic 
distribution is characterized by the appearence of the 
tilting mode that emittances are periodically exchanged 
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FIG. 2: Evolution of the emittance (top) obtained of the 
equation ((5]). ex and ey are represented by red and blue 
lines, respectively. Emittance and sum transversal emittance 
(bottom) as a function of the ratio of the envelope beam 
X~ f^l'^'v ^x, ey and (e^, + ej,)/2 versus x s^re represented by 
red, blue and gray lines, respectively. 



between x and y, similar to a second order difference res- 
onance driven by skew quadrupoles. We note that this 
second order difference resonance is entirely absent, if the 
initial emittances are chosen equal. This tilting insta- 
bility between x and y obviously re quir es a sufSciently 
anisotropy. Note that according to (22| the symmetric 
perturbation only excites the breathing mode, whereas 
the antisymmetric perturbation excites the quadrupole 
mode only. With anisotropy this is not the case and we 
find mixing of breathing and quadrupole modes. In our 
model both modes are resonance. Space-charge induces a 
coherent shift of the resonance conditions luximujy = 0, 
since the full ensemble of particles respond to the reso- 
nance in a coherent way. For such a coherently oscillat- 
ing beam along an uniform linear focusing channel the 
resonance condition is lujx ± mujy + Aoj = p^ . For 
the tilting mode we expect different shifts for the cases 
ujx+ujy-\- Alo = (sum resonance) and — tjj, -1- Auj = 
(difference resonance) [l^l . Auj is the coherent shift away 
from the single particle resonance condition caused by 



the coherent motion of all, or a large fraction of particles 
of beam. It makes the majority of 2500 test particles 
launched along the x- and y-axes inside of the beam to 
be resonance as shown by the histogram in Fig. 3) In our 
model the driving term for this difference resonance is not 
a skew quadrupole as in synchrotrons, but the internal 
space-charge force caused by the exponentially growing 
tilting of the beam cross section. There is a strong de- 
pendence of the difference resonance on the emittance 
ratio. In cases where the difference ressonance is crossed 
for unequal emittance we expect an asymmetric behavior 
of beam. 
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FIG. 3: Asymtotic evolution of the envelopes (top) and emit- 
tance (bottom) obtained of the equations (|3| and ©, respec- 
tively, (top), ex (bottom) and ry (top) ey (bottom) are 
represented by red and blue lines, respectively. 



Nonlinear resonances eventually yield rms emittance 
growth as more and more particles to be launched out of 
core. Therefore, the halo is not caused by a collective ef- 
fect involving all core particles but is caused by resonant 
interaction of some particles with the mismatch modes. 
Assuming which beam is usually equipartitioned in x and 
2/ = {''"y^x)'^ /{rx^yy), where ^ is the ratio of oscilla- 
tions energies in the x and y directions [s^l, but it has 
large beam size-rms mismatched, the resonances enable 
"excess" energy transfer from one plane to another p^ . 
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We show that the exchange is accompanied by halo cre- 
ation along one direction preferential. 



FIG. 4: Histogram of the dimensionless frequencies uix and Uy 
of 2500 test particles. We launch the test particles along the 
X- and j/-axes of beam in specified region inside the beam (be- 
tween O.Olrj; and Q.Itx spaced by 0.0004ra; along a;-direction 
and between O.Olry and 0.7rj; spaced by 0.0004r-j; along y- 
direction). These 2500 test particles has Ux = 1.05461 and 
2455 test particles has ujy — 1.05461 and 45 test particles has 
iUy = 1.00667. We used FFTW (The Fastest Fourier Trans- 
form in the West) [58l | to compute numerically the frequencies. 

To analyse the resonances we computed dimensionless 
frequencies of breathing and quadrupole modes and the 
dimensionless frequencies cOx and iOy of 2500 test particles 
launched along the x- and y-axes inside of the beam using 
Fourier analyses. We used FFTW (The Fastest Fourier 
Transform in the West) to compute numerically the fre- 
quencies. Breathing mode frequency is ujx^ = 2.01334 
and quadrupole mode frequency is ujx^ = 2.01334. Note 
that both modes has the same frequency. Therefore 
breathing mode and quadrupole mode are resonance 
(1 : 1). 2500 test particles has lu^ = 1.05461 and 2455 
test particles has ujy — 1.05461. Thus 2455 test parti- 
cles obey resonance condition oj^ — i^y = 0, and 45 test 
particles has LUy = 1.00667 as shown by the histogram 
in Fig. [31 The system is dominated by the ressonance 
ujx — uJy = 0, which creates a barrier between the re gion 
inside the beam and the region outside the beam [49|. 
The space-charge coupling force is responsible for the 
energy transfer from the core oscillations to the single 
particle oscillations. The dominant \th order ressonance 
at u>x/u>y = 1 yields significant emittance exchange [4^ . 
The particle-core (2:1) resonances are present in both 
quadrupole mode and breathing mode of core oscillation, 
but only in y-direction of test particle [2l|. Thus, halo is 
formed along y-direction as illustrated in the Fig. [5] [23| . 
The key to interpreting anisotropic halo growth due to 
the (2 : 1) is the dependence of space-charge force, which 
increases as the beam size becomes smaller. The para- 
metric (2 : 1) resonance halo, modified by driving term 
for the difference resonance that is the internal space- 
charge force caused by the exponentially growing tilting 
of the beam cross section, remains the dominat mecha- 
nism to explain halo formation along one direction pref- 
erential. 

For large beam size-rms mismatched and initial ra- 
tio envelopes beam x ~ li the ratio of oscillations 
energies in the x and y directions remaines constant, 
5 = 1 [S^]. As illustrated in the Fig. [5] beam fast suf- 




FIG. 5: Evolution of the dynamical test particle, x — Q.Sr^ 
and y — O.bVy, initial. The orbits x and y of test particle are 
represented by black and pink lines, respectively. The beam 
envelopes, Vx and ry are represented by red and blue lines, 
respectively. 



fers anisotropization characterized for discontinuous vari- 
ations in X s-nd r] |23|. l56| . The anisotropy leading to 
coupling resonance [171 . Il8l | in the presence of nonlinear 
space-charge forces was suggested as a possible approach 
to the equipartitioning question (26l . [27l [3ll [3a |. since 
collisions cannot be made responsible for energy transfer 
in linacs. The underlying mechanism is collective oscilla- 
tions of the space-charge density that creates nonlinear 
forces similar to those magnetic sextupoles and leads to 
the resonat coupling [13, H^, [l^l ■ 




FIG. 6: The ratio of oscillations energies ^ and anisotropic 
ratios, x ^'Hd rj in beam propagation. ^ is represented by red 
line, and x ^^d rj are overlapping represented by blue line. 

Statistical theories of chaotic systems are generally 
based on Liouville's theorem and the assumption of 
eqnipartition, subject to the constraints given by invari- 
ants of the system. If the only invariants of the system 
are the total energy and the total number of particles, the 
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assumption of equipartition leads to the theory of statisti- 
cal mechanics at thermal equilibrium. The temperature 
is then uniform. In turbulence, the energy is not con- 
served, since there is an external energy source, but there 
may be other invariants that survice the turbulence. The 
random flutuations then drive the system toward turbu- 
lent equipartition (TEP) on the hypersurface in phase 
space that is defined by these invariants (85j . Previous 
examples of TEP in plasmas physics were considered by 
Yankov [s^l and Isichenko et al. [s^ . Turbulent equipar- 
tition corresponds to a phase space density uniform on 
the surfaces of the two constant adiabatic invariants, a 
version of the ergodic hypothesis where adiabatics invari- 
ants play the role of the no longer conserved energy. The 
applicability of the statistical mechanics to self-consistent 
plasma turbulence is not clear, but some predictions of 
this approach are in a qualitative agreement with exper- 
iments 'S^. In the state of turbulent equipartition, the 
plasma density and temperature are inhomogeneous, and 
there are no particle or energy fluxes. Such a state is usu- 
ally marginally stable to the modes potentially responsi- 
ble for driving the plasma to this state. We believe that 
the problem of anisotropic equipartition in linac should 
be solved in the same spirit. 



III. ANISOTROPIC KINETIC AND DYNAMIC 
PROCESSES IN EQUIPARTITIONED BEAMS 

In order to understand the initial dynamical behavior 
of an anisotropic beams, in particular, to study possi- 
ble mechanisms of equipartition connected with phase 
space we have to know how we can compute the vari- 
ables (volume, area of surface, and area projected) that 
characterize the anisotropic beam in phase space. The 
equipartition is generally based on the assumption that 
system is ergodic. If the dynamics is chaotic in some 
subspaces of the phase space, the system can be consid- 
ered as ergodic in these subspaces. The beam began in a 
nonequilibrium state evolved toward a metaequilibrium 
in which the particle orbits filled an invariant measure of 
phase space. The transient dynamics reflected an intri- 
cated network of space-charge modes. We shall call this 
subspace of the phase space " C " i where ^ is the ratio of 
oscillations energies in the x and y directions. Unlike the 
thermodynamic equipartition energy in conservative sys- 
tems, anisotropic equipartitions describe strongly non- 
equilibrium systems such as high-intensity ions beam in 
space-charge dominated regime with large mismatched 
beam size-rms initial. Anisotropic equipartition corre- 
sponds to a phase space density uniform on the surface 
invariant of the — [ryexY / {rx^y)"^), ^ version of the er- 
godic hypothesis where the ^ invariant play the role of the 
conserved energy [g^] ■ In the state of anisotropic equipar- 
tition, the Lyapunov function [o^ and temperature are 
stationaries, the entropy grows in the cascade form, there 
is a coupling of transversal emittance, the beam develops 
an elliptical shape with a increase in its size along one 



direction and there is halo formation along one direction 
preferential. Such a state is usually marginally stable. 

Emittance is the area in phase space occupied by the 
particles of the beam. The maximum emittance of a 
beam that a system can accept is called the acceptance 
of the system " a " [sO, [Mi ■ Then the acceptance is 
defined as the maximum phase space area where particles 
can survive in the accelerator. It can be determined by 
dynamical effects associated to space-charge forces. In 
a linear accelerator with nonlinearities, the transversals 
motions are coupled. The particles move on the distorced 
surface in four-dimensional transverse phase space, the 
so-called "hyper-egg" . We can talk only about the pro- 
jections of this hyper-egg onto the two transverse phase 
planes ax{x,x ) and ay{y,y ). These projections are no 
longer clean curves but bands. It is well-know fact that in 
presence of space-charge the beam phase space structure 
gets quite complicated and is divided to areas of stable 
and unstable motion. 

Qualitatively, emittance can be thought of as the tem- 
perature of the beam, a measure of the random disor- 
der in the transverse momenta of the constituent parti- 
cles. Quantitatively, transverse emittance is defined by 
drawing a contour around a given percentage of particles 
in phase space and gives us a numerical figure of merit 
for describing the quality of the beam. As the beam is 
focused or defocused, the convergence or divergence of 
the beam envelope yields a correlation between particle 
position and angle of motion; if, however, we focus the 
beam to a waist, the correlation between particle posi- 
tion and transverse momentum is minimized. We de- 
fine the RMS acceptance as the maximum square root of 
the ttj — — {jj')^ (henceforth, j ranges over 

both X and y, where x and y are the positions of the 
beam particles). The second moment term {jj )^ repre- 
sents a correlation between j and j that exists when the 
beam envelope is converging or diverging; qualitatively, 
it can be thought of as a measure of inward or outward 
flow of transverse kinetic energy. At a waist, this cor- 
relation is minimized and the second moment term is 
zero [53!]. Thus, the acceptance reduces to maximum of 
the Qj = \/{p){f^- It is a very important for one to 
know how much of the phase space is stable. With the 
help of the Lyapunov functions we can construct an in- 
variant beam area. This function is dependent of the 
equipartition ^, of the anisotropy variable 77, ratio of the 
average external focusing force to the space-charge force 
a and space-charge strength p. 

It is intuitively clear that if the total energy of a phys- 
ical system has a local minimum at a certain equilib- 
rium point, then that point is stable. This idea was 
generalized by Lyapunov [93] into a simple but power- 
ful method for studying stability problems in a broader 
context. Lyapunov functions are functions which can be 
used to prove the stability of a certain fixed point in a dy- 
namical system or autonomous differential equation. For 
dynamical systems (e.g. physical systems), conservation 
laws can often be used to construct a Lyapunov func- 
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tion. A certain subclass of dynamical systems, namely 
potential or gradient system are of particular interest be- 
cause their behaviour is simpler than the general case, 
and because they are frequently encountered in approx- 
imate treatments of physical systems. For a gradient 
system, = f{x), takes the form = 

where L{x) is a function of the variable x. More gener- 
ally, a system with an equilibrium xq is said to have a 
Lyapunov function L for this equilibrium if this function 



satisfies the conditions L{xq) = 0, 



dL{x) 
dt 



< for X ^ Xq 



and L is a smooth function of x in some neighborhood of 
Xq [HI] . A gradient system, satisfying ^^^^ 



dL(x) 
dx 



has 



global Lyapunov function if L is bounded. For gradient 
system the dynamic consists of relaxation toward mini- 
mum in L. This means that such functions are strictly 
only defined when corresponding equilibriums are fixed 
points. Finding a Lyapunov function for a certain equi- 
librium might be a matter of luck. Trial and error is the 
method to apply, when testing Lyapunov functions on 
some equilibrium. 

Next, we will apply Lyapunov's method above to the 
acceptance dynamics [9J]. To this end, we will replace 
the aforementioned quantities x(t) and L(x) by the aj{s) 
acceptance, and L(dj) Lyapunov function of the accep- 
tance. A definition of a Lyapunov function would be : 



L 



a 1 




p{l + l/v) 



(6) 



where aj (henceforth, j ranges over both x and y, where 
X and y are the positions of the beam particles) is accep- 
tance in the absence of directional correlations between 
j and j , primes denote derivates with respect to s. This 
function is dependent of the beam envelope rj, of the 
equipartition ^, of the anisotropy variable 77, ratio of the 
average external focusing force to the space-charge force 



a and space-charge strength p. To a^g — and ay 



0, 



L satisfies the condition L(axo,ayg) = 0. It can be in- 
terpreted as the equation of a surface that resembles a 
paraboloid opening downward and tangent to the Qxtty 
plane at the origin. L can be considered as a "free en- 
ergy " [q^I ■ To the issue '^''"'j"^ < mentioned above the 
evolution of the L{dj) Lyapunov function represented by 
blue line (top graphic) in the Fig. [5] can be regarded as 
a proof. As the Lyapunov function ([6]) satisfies the con- 
ditions L(dj„) — and '^"^j"^ < we can apply equation 
"~ds~^ = ——^^ obtaining the following acceptance dy- 
namics equations: 



ddx 
ds 
da 



^ = -a i-yT7eP^ + p(i + i/^)5. (7) 



where s is the axial coordinate of a beam. Equations ([7]) 
and ([5]) express the acceptance dynamics in term of 
the beam envelope r^, of the equipartition ^, of the 
anisotropy variable 77, ratio of the average external fo- 
cusing force to the space-charge force a and space-charge 
strength p. We derive equations for acceptance change 
in each plane for continuous elliptical beam within a con- 
stant focusing channel without correlation between par- 
ticle position and transverse momentum. The equation 
will apply to beams within a conducting pipe whose ra- 
dius is much larger than the beam size. The resulting 
acceptance equations contain two terms: the first term 
describe acceptance changes associated with transfer of 
energy between the two planes; the second describes ac- 
ceptance changes associated with anisotropic processes. 
The anisotropy leads to coupling resonance [13, [H, [111 
in the presence of nonlinear space-charge forces. Space- 
charge couple some degrees of freedom causing resonance 
between them. Coupling resonance leading to an ex- 
change of energy between the degrees of freedom. 

The general solutions to the equations ([7]) and ^ are 
given by: 



cie 



C2e 



■f (\/i7«-i) +p(i+i/-'()] ds 

l[T^{V^-l)+p{l+ri)]ds 



(9) 
(10) 



where ci and C2 are arbitrary constants determined by 
initial values of the solutions ([9]) and (fTO]) . respectively. 

The oscillations beam envelope perturb nonlinear 
space-charge force yields a correlation between parti- 
cle position and transverse momentum. As the beam 
compresses and expands the particles gaining and losing 
kinetic energy. Thus, the second moment term {jj )^ 
is not zero and the RMS acceptance becomes aj = 
'\/0'^)0''^) ~ {jj')'^- Therefore the equipartition and the 
variable anisotropy are given by ^ = {ryaxY / {rxCLyf' ^^'^ 
V = cLx/cLy, respectively. Thus the equations ([7]) and ([8]) 
are transformed by : 



dax _ 
ds 
day 



+ p] ax + { ^ + p] ay, (11) 



— +p]ay (12) 



where the terms axjvx and ayjvy model the equiparti- 
toning and the term ax + ay covers residual growth from 
nonlinear resonance 38j. Jameson has derived identical 
equations [H, [H, [34 1 to the question of equipartioning 
in linear accelerator. 

To solve the equations (HH) and (O we add (fTTI) 



and and then we solve for ax^ay — cae'^''*'. Then 
it is easiest to substract ((T^ from (fTTI) and substitute 
to obtain the general solutions to 



^~ \ a ^ ^"^ ~ ~'^y ^ cse'-^'"'-' to obtain the 

-T^ a \\fl - 1 ) — + P (1 + ??) dy (8) the equations ((n|) and given by 

GjS \ / T'u 
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ax = cse^'"' - \cz 



ay = < C3 



ds + C4 > < e 



(is + C4 > < e 



(13) 



(14) 



where C3 and C4 are arbitrary constants determined by 
initial values of the solutions ([13]) and (|14p . respectively, 
aj is dependent of the beam envelope dynamic 7'j, of 
the anisotropy variable variations x = ^x/fy and space- 
charge effects, a and p. 

To analise the effect of coupling between the two trans- 
verse phase planes without {ax{x,x ) and ay{y,y )), and 
with {ax{x,x ) and ay{y,y )) correlation between parti- 
cle position and transverse momentum (the second mo- 
ment term {jj )^ = and {jj )^ ^ 0, respectively) we 
solve the integrals in © and ^ and HH) to ob- 

tain the acceptance evolution without and with correla- 
tion, respectively. We launch the beam with a — 0.00001 
and p = 0.25 in space-charge dominated regime, mean- 
ing that the collective oscillations dominate over the in- 
dividual particles betatron motion, ci — 0.22360 and 
C2 = 0.22360, C3 = 0.44721 and C4 = 0.22360 are de- 
termined by initial values dxo — o,yQ — 0.22360, and 
a-xo = o,yo — 0.22360 of the solutions ([9]) and pO)) . 
and and (fH|) . respectively. We solve the integrals 
in (HI) and (Uni), 1131) and dUl) up to s = 50.0. The cor- 
responding evolution of the rms-emittance ([S]) and the 
rms-envelope ([4]) are used. 
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FIG. 7: Evolution of the beam acceptance without [hx and 
5y), and with (a^ and ay) correlation between particle posi- 
tion and transverse momentum, and ay are represented by 
pink and black lines, respectively, a^ and ay are represented 
by red and blue lines, respectively. 



The acceptance evolutions without and with correla- 
tion between particle position and transverse momen- 
tum given by equation ([9]) and (fTO|) . ([TS]) and ([T4|) . re- 
spectively are shown in Fig. [T] The area in a; — x and 
y—y trace spaces defined by hx (pink lines) and hy (black 
lines), respectively, remain constant with very small os- 
cillations because hx and hy are beam area projections 
derived from Lyapunov functions ([5]). Also, if a energy- 
balanced beam is injected, the distribution will move to- 
ward equipartition, in which the beam assumes a uniform 
distribution on the surface of cix and ciy [s^ . However, it 
is observed transverse phase planes coupling in a; — a; and 
y — y trace spaces defined by ax (red lines) and ay (blue 
lines) caused for space charge driven core-core resonance 
together with single-particle resonances (T^. [30l[38|.This 
coupling is characterized by the energy exchange be- 
tween the directions. Energy/acceptance exchange re- 
quires resonant coupling, which can take place only if 
an intrinsic resonance relationship is fulfilled. A signifi- 
cant number of particles are trapped inside a resonance 
island [49, ^] as shown in Fig H) They are characteriz- 
ing the formation of a metastable state of the particles 
beam (s^, [11] . Geometrical methods can be used to an- 
alyze this state [H, [S^ . By definition of acceptance we 
observe the similarity with the evolution of emittance {cx 
and Cy) shown in Fig.[2l 

In the process of acceptance transfer, strong correla- 
tions between particle position and transverse momen- 
tum are naturally developed. The acceptance oscillations 
are now driven by variations of the space-charge force due 
to the beam compresses and expands. The rapid acep- 
tance oscillations are due to coherent transverse plasma 
oscillations in the beam and are a manifestation of peri- 
odic energy exchange between potential and kinetic en- 
ergies. The acceptance growth results from nonlineari- 
ties in the particles oscillations about their equilibrium. 
Excess energy is required for driving the nonlinearities, 
which stems from the energy anisotropy between differ- 
ent degrees of freedom. Acceptance grows in a plane 
that receives energy, and decreases in a plane that loses 
energy. The acceptance ax growth increased, and the ac- 
ceptance ay growth decreased and vice-versa during the 
propagation of the beam. Such transfers between x — x 
and y — y are commonly observed if the initial accep- 
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tance differ. This indicates that ratios acceptance are 
important variables. In order to expedite an acceptance 
flow from one direction to another, the two directions 
have to be on couphng resonance. TypicaUy, we employ 
difference resonance: luJx — mtOy — where lu^ is single 
particle frequencies along of a;-direction and ujy is single 
particle frequencies along of y-direction, I and m are inte- 
ger numbers. In our case I = m = 1 and number large of 
particles obey resonance condition lOx ~ "-^y = shown 
in Fig.m Acceptances only start to flow from a "hot "to 
a "cold "degree of freedom. This suggests the develop- 
ment of correlation among the degrees of freedom during 
the acceptance exchanging process. 

Energy/acceptance flow is due to the correlation jj 
between the velocity and the position of the particle in 
regions where the beam size contracts or expands. The 
term {jj )^ represents an inward or outward flow term 
in the transverse kinetic energy. The physical interpreta- 
tion is that the rms kinetic enrgy of the particle distribu- 
tion consists of a thermal component and a flow compo- 
nent [s^l ■ Transfer of kinetic energy from one coordinate 
direction to another results in partial or complete kinetic 
energy equipartitioning. While superficially similar to 
equipartitioning of energy in a gas, this transfer does not 



result from individual particle collisions but, presumably, 
from interactions between the individual particles and 
collective fields. Hofmann has show that coherent-mode 
instabilities can lead to kinetic-energy exchange 



17| 



In subsequent analysis, the ratio a^jay is varied to ex- 
plore different degrees of anisotropy and energy flow. The 
rate of exchange can be predicted with certainty analyt- 
ically. To obtain the relation between a^, a-nd 5,^, and 
and Oy, respectively, which will give the acceptance flow 
from one direction to another, elimanate s by dividing ([7]) 
by (HI), and (fTTjl by (fT2|) . respectively. The equations for 
the acceptance flow are: 









dx 


ddx 




ddy 






dy 









(15) 



(16) 



The equations (fTo)) and p6|) are separable variables. 
The general solutions are : 



dax {—a + prx)ryax + (a + pry) 
day (a -I- prx) fyax -\- { a -\- pfy) Txay 



F- i~ 



F, 



Oiy j Oi'j; 



ryln {ax) - In [ay] , 
6i In (bge^ + 646 + 65) / 2^2 
26^ ^ 



b3 



arctan 



2636 + bj 



- In (ax) 



(17) 
(18) 



where 77 = e^/ey, 6 = Ux/ay, bi = -ry {pvx + a), 
b2 = Tx {pry - a), 63 = ry {prx + a), 64 = a {rx - ry), 
bz — rx {pry + a). rx,ry is the beam envelope, a is the 
ratio of the average external focusing force to the space- 
charge force and p is the space-charge strength. Fg^/g^ 
is the acceptance flow without correlation and Fa^/a^ is 
the acceptance flow with correlation. 

The necessary conditions for efficient heat transfer be- 
tween transverse degrees of freedom are formulated. As 
discussed, {jj )^ represents an inward or outward flow 
term in the tranverse kinetic energy. The tranverse ki- 
netic energy consists of a thermal component and a flow 
component. The latter is due to the correlation between 
the velocity and the position of the particle in regions 
where the beam size contracts or expands. Flow com- 
ponent is then directly related to heat transfers between 
different degrees of freedom within the beam. The dif- 
ference between the acceptance flow with correlation and 
the acceptance flow without correlation is the flow com- 
ponent defined by F^^' = F^^/a, - -f^a^/a, ■ The flow com- 
ponent evolution F-^i is shown in Fig. [Sj It is affected 
by a swift variation similar discontinuous variations in x 
and 77 in Fig [5] The amplitude oscillations grow during 



the propagation of the beam. When the beam radius ex- 
pands, the flow component has an outward direction and 
the thermal component decreases, when the beam radius 
contracts, the flow component is inward and the thermal 
component increases. The propagation of the beam is 
thus characterized by a variation in the envelope beam 
which correlates with the generation of the heat flow and 
a variation of the beam temperature. As usual in sta- 
tistical physics, we relate the temperature of a particle 
ensemble to its random motion. Charged particle beams 
change their size while passing through an ion optical 
system. The tranverse kinetic energy (p|)/2m contains 

a flow component if {jj ) ^ 0. Therefore, the flow com- 
ponent of the kinetic energy must be subtracted from 
transverse kinetic energy in order to obtain its thermal 
component. The "non-equilibrium temperature " Tj is 
defined as the thermal component of the transverse ki- 
netic energy of the j-th degree of freedom. 

The elementary mechanism responsible for the transfer 
of heat from one degree of freedom to another is consti- 
tuted by the effect of interactions between the individual 
particles and waves of collective fields [5^. The trans- 
fer of heat arises from the scattering of a particles and 
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Evolution of the flow component F. .' obtained of 



the difference between equation (|18|l and the equation (I17p . 
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waves by a mean-background of waves of collective fields 
in thermal equilibrium via resonant interactions [t^ . [t^ . 
We demonstrate that resonant interactions of the parti- 
cles beams provide a mechanism for effective energy ex- 
changes among different degree of freedom. The issue of 
reaching thermal equilibrium is to imagine a heat bath, 
which is in thermal contact with the beam. In this case, 
the particles of the beam can exchange energy through 
the thermal bath. Chavanis [6J] uses a kinetic theory to 
describe the relaxation of a test particle in a thermal bath 
of field particles. The relaxation of the system is due to a 
condition of resonance and it may happen that the relax- 
ation stops because there is no resonance anymorc.lt has 
been shown that all heat transfers within the beam feed- 
ing thermal energy from one degree of freedom to another 
are always associated with an increase of beam entropy 
and thus always lead to an irreversible degradation of 
beam quality. Since the particle distribution of a beam 
is confined by focusing potentials and the individual par- 
ticles are performing oscillations, there is a continuous 
exchange between potential energy and kinetic energy so 
that displacements in position due to random processes 
translate into velocity changes and vice-versa. The ba- 
sic idea is that a continuous exchange between potential 
energy and kinetic energy can trigger chaos via a para- 
metric resonance [4, 60]. 

A priori there need be no direct connection between 
increases in chaos and exchange of energies. However, 
one would expect resonant couplings to lead significant 
changes in energies. A correlation between changes in 
energy and the amount of chaos thus corroborates the 
interpretation that this chaos is resonant in origin. This 
transient chaos can drive chaotic phase mixing, which, 
in the context of a fully self-consistent evolution, might 
be expected to play an important role in violent relax- 
ation [gl, [63, HI] ■ After the energy redistributes among 
all the modes to achieve "thermal" equilibration [6^1 ■ Re- 



cent simulations (28| have shown that fully self-consistent 
simulations of beams can exhibit evidence of chaotic 
phase mixing. 

To characterize the "thermal "equilibrium we analyze 
the dynamics of the Lyapunov function ([6]), and of the 
temperature and entropy of the beam. The "temperature 
"Tj of the j-th degree of freedom of a charged particle 
beam can be expressed in terms of second order beam 

moments T, ~ t—t-, where m is the mass ions beam 

and kf, is the Boltzmann's constant [gl, [T^I- This for- 
mula shows that envelope Vj and acceptance variations 
cause temperature variations. For a system with "non- 
equilibrium temperatures "that is oscillating anisotrop- 
icly around Teg, the equilibrium temperature can there- 
fore be aproximated by the arithmetric average of the 



T. 



eq 




(19) 



With the equilibrium temperature Teg for the 2-D 
beam model, the entropy change near thermodynamic 
equilibrium due to a temperature balancing process may 
then be written as 



ds 



TxTy 



(20) 



Obviously, the entropy S{s) remains unchanged in the 
case of temperature equilibrium while increasing during 
temperature balancing. We may regard Eq.(l20|) as a 
particular manifestation of Boltzmann's H-theorem [7l| . 
The basis for the dynamics behaviour of the entropy is 
the relation between the acceptance, the beam tempera- 
ture and the envelope. Within the beam, heat exchange 
between the degrees of freedom may occur, leading to 
an entropy growth as described by Eg. ([^0)1 . We con- 
clude that equipartitioning effects occurring within ini- 
tially thermally unbalanced charged particle beams are 
always associated with an irreversible degradation of the 
beam quality as a whole. Beam transport without an 
increase of entropy are thus possible if either the beam 
stays round throughout its propagation. The entropy 
concept was first applied to beam, and its relationship to 
rms emittance was explored by Lawson, LapostoUe and 
Gluckstern [73 |. In 1991, the maximum-entropy hypoth- 
esis was used to calculate the characteristic of the final 
distribution for a high-intensity expanding beam in free 

space Its']. 

We launch the beam with m — 1, rjo — 0.43616 and 
fljo — 0, 22360, consequently Tjo = 0.26282 initially, and 
integrate the entropy equation ([20]) up to s = 200.0. The 
corresponding evolution of the temperature Tj is used 
in Eq. (|^D|) . To the evolution of the temperature Tj the 
corresponding evolution of the rms-acceptance (jlip and 



([T2)l . and the rms-envelope Q are used. To the Lya- 
punov function dynamics the corresponding evolution of 
the rms-acceptance ([7]) and dH) and the rms-envelope ^ 
are used. The dynamics of the entropy S (bottom), tem- 
perature Teq (top) and Lyapunov function L (top) are 
shown in Fig. [9l 
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FIG. 9: Evolutions of the temperature T^q (top), Lyapunov 
function L (top) and entropy S (bottom) obtained of the 
equations (IT^ . © and (PP)) . respectively. Teq (top) and 
S (bottom) are represented by red lines and L (top) is repre- 
sented by blue line. 

In the case without correlation between particle posi- 
tion and transverse momentum the Lyapunov function L 
is a monotonically decreasing function with respect to s 
for solutions Oj (because we have '^^Ij^' < 0) as shown in 
Fig. [9l the fixed points cixo = and dyg = correspond 
to extrema of L (i.e. we have = '^^Jh' ^ ~ 

and as L is bounded then from the two aforementioned 
properties it follows that L becomes stationary in the 
limit s ^ oo. The stationarity of the Lyapunov function 
L implies the stationarity of the ions beam. 

In the case with correlation between particle position 
and transverse momentum the temperature T^q oscillates 
after steady with small fluctuations and entropy S grows 
by cascades as shown in Fig. (5] By analogy with com- 
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pression and expansion of a gas, the beam temperatute 
Teq heats up during compression and cools during expan- 
sion. The concept of entropy cascade is the key agent in 
the heating and relaxation of the beam ^TG'l.The phys- 
ical process by which this relaxation occur is the wave- 
particle interaction rather than the familiar Coulomb col- 
lisions: the reservoir of free energy, represented by the 
Tx/Ty > 1 or Tx/Ty < 1 states, pump energy into ini- 
tially small microscopic fluctuations in the system and 
lead to the build up of a signiflcant level of electromag- 
netic anisotropic flelds. This is done initially in the tem- 
perature Teq which therefore decreases in time. As a 
non-linear feedback effect, the waves heat the particles 
of the system. The feedback effect of the waves on the 
particles results in a final state in which particles and 
fields are in a quasi-stationary state [t^, uM- Qualita- 
tively, the coUisionless relaxation is driven by the fluc- 
tuations of the fleld, itself induced by the fluctuations 
of distribution function. The fluctuations of the electro- 
magnetic flelds are able to redistribute energy between 
particles beam and provide an effective relaxation mech- 
anism on a very short timescale. Physically, this coUi- 
sionless relaxation is interpreted as reflecting a resonant 
coupling of particles with the wave of flelds The 
resulting "resonant phase mixing "might be sufficiently 
strong to explain violent relaxation. During violent re- 
laxation, the beam tends to maximize the rate of entropy 
production ^ while conserving the constraints imposed 
by the dynamics. Moreover, as the system approaches 
quasi-equilibrium, the fluctuations of the fleld are less 
and less efficient. 

CoUisionless Landau damping of the electromagnetic 
fluctuations leads to particle heating in the sense that 
it transfers the electromagnetic fluctuation energy into 
fluctuations of the particle distribution function, which 
are then converted into heat by stochastic space-charge 
forces 8Qj . The nonlinear collective forces have the same 
effect as collisions in thermalizing a particle distribution. 
The available energy is, however, not entirely thermal- 
ized. Some of the energy will be converted to poten- 
tial energy, due to the change in beam radius. Stochas- 
tic space-charge forces are required to increase the en- 
tropy. The entropy cascade is the way in which the 
energy diverted from the electromagnetic fluctuations 
by the coUisionless damping (wave-particle interaction) 
can be transferred to the stochastic space-charge forces 
scales resulting in equipartitioning. Equipartitioning via 
this anisotropic instabilities can properly be described in 
terms of a thermalization of the beam in the sense of dif- 
fusion riii, The beam distribution never reaches 
the true equilibrium distribution. Instead, two distinct 
distributions are seen the majority of the beam is very 
nearly in equilibrium, surrounded by a smaller density 
beam halo. We deflne the halo as the nonthermal part of 
the beam that consists of particles having oscillation am- 
plitudes that are larger than the maximum of the beam 
core. 

We observe to have different regimes depending on the 
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value of time scales. There is first a phase of violent relax- 
ation on a time scale s = 50 leading to a quasi-stationary 
state. This phase is followed by a thermalization leading 
to the Teq stationary on a time scale s = 120 due to the 
thermal bath (field waves) i.e^^he combined effect of 
imposed friction and diffusion [sil, IH, [s^ , the diffusion 
process are arising from the fluctuations of the self-fields 
and the friction results from a polarization process; not 
to collisions effects. The first phase is described by the 
Vlasov-Poisson system with building blocks of the beam 
core in phase space equipartitioning [gl, [s^l and the sec- 
ond phase by the Lenard-Balescu-Poisson system [82, 's^. 
But there is a time-scale separation between the phase of 
violent relaxation and the phase of thermal bath relax- 
ation. We can consider for intermediate times s = 50 
to s = 120 that the distribution function is a quasis- 
tationary solution of the Vlasov equation of the form 
/ = f{E,t) that slowly evolves under the action of im- 
posed friction and diffusion thermal bath, not collisions. 
Therefore, the system first reaches a state of mechanical 
equilibrium through violent relaxation , then a state of 
thermal equilibrium through the effect of imposed fluctu- 
ation and dissipation i.e., the thermal bath. The study 
this dynamic will be considered in a future work. 

The dynamics of particles described by these processes 
has a complex phase space structure. The effectively ac- 
cessible phase space of the beam can have a complicated 
geometrical structure. Qualitatively, the accessible phase 
space of the beam can be computed by volume and sur- 
face area that the beam occupies in phase space. We 
can then measure the volume and surface area of phase 
space spanned by beam. A commonly employed measure 
of volume is the root-mean-square geometric emittance 
and of the surface area is the sum of the 



emittance a — ex + Sy [89|. To the parabolic beam of 
the Sec. the volume ^ and surface area a are calcu- 
lated analytically using the emittance equation ([5]). The 
surface area a is given by : 



cr = —VT5KA 
90 



(x' + 1) (Sx' + 2x + 5) 



(i + xr 



1/2 



(21) 



and the volume fi is given by 



n = —VT5KA 
^ 90 



(5x^ + 2x^ + 5x) 

(1 + x)' 



-I 1/2 



(22) 



where K is the dimensionless perveance of the beam, A 
is mismatched amplitude of the oscillatories modes beam 
and X = '''x/ry is ratio of the envelope beam. The tem- 
poral evolution of the volume /i and surface area cr of the 
beam are shown in Fig. [TOl 

As shown in Fig. [TU] the dynamics of fi and a present 
jumps. They are coupled because when the evolution of 
the volume /x is maximum, the evolution of the surface 
area ?7 is a minimum and vice versa, -i.e. the particles 
oscillate between the center and surface of the beam in 
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FIG. 10: Evolution of the volume and surface area obtained 
of the equations (|2H) and (|22[) . respectively, fi and a are 
represented by red and blue lines, respectively. 



real space X — Y and in the phase space, when the evo- 
lution of the surface area is minimal, the particles are 
distributed along the volume and when the evolution of 
the surface area is maximum, the particles are distributed 
on the surface gaining velocity. 

In the Fig. (jlip the volume fj, and surface area a vary 
in function of x- For x < 1 A* increases and a decreases. 
At X = 1 the surface area has the minimum value and the 
volume has maximum value. For x > f M decreases and 
a increases. At x ~ 150 surface area and volume become 
stationaries . Note that the surface area is always greater 
than the volume in function of anisotropy variable x- 
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FIG. 11: Volume and surface area as a function of anisotropy 
variable x = '''x/ry. a and n, versus x ^-re represented by red 
and blue lines, respectively. 

The volume /i and surface area a are not constant in 
function of the axial coordinate of the beam s. But 
the ratio of oscillations energies in the x and y direc- 
tions ^ remaines constant in function of the s as illus- 
trated in the Fig. [51 The ^-space is the relevant one 
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to make the statistical mechanics of anisotropic beam. 
Arguably, the only crucial point is the possibility of con- 
structing time-independent solutions to the Vlavov equa- 
tion which depend on quantities other than global iso- 
lating integrals such as energy, so as to ensure that, if 
initial data be evolved into the future along the charac- 
teriztics associated with the self-consistent potential, the 
form of initial distribution function /o will remain un- 
changed [sljIU.The key point is that, in principle, a self- 
consistent equilibrium can be constructed from any set 
of time-independent phase-space building blocks which, 
when combined, generate the density distribution associ- 
ated with an assumed time-independent potential. 

As the surface area and volume become stationaries 
in function of the anisotropy variable x we can develop 
an approach to modelling locally anisotropic kinetic pro- 
cesses [131 ■ The generalization of kinetic equations to 
anisotropics bams is not a trivial task. In order to con- 
sider the anisotropic processes of a microscopic or macro- 
scopic nature it is necessary to reformulate the kinetic 
theory in anholonomic frames. The modelling of ki- 
netic processes with respect to anholonomic frames is 
very useful with the aim of elucidating flows of fluids 
of particles not being in local equilibrium. The deduc- 
tion of the equation for one-particle distribution function 
in anisotropics phase-space will be considered in a future 
work. 



IV. CONCLUSIONS 

We examined the transition from isotropic to 
anisotropic beam profiles in a linear focusing channel. 
We considered a high-intensity ion beam in space-charge 
dominated regime and large mismatched beam size-rms 
initial and we observed a fast anisotropy situation of 
the beam characterized for a transition by the transver- 
sal section from round to elliptical with the coupling of 
transversal emittance. 

The coupled motion between the two tranverse coordi- 
nates of a particles beam is arising from the space-charge 
forces. A beam with nonuniform charge distribution al- 
ways gives rise to coupled motion. However, it is only 
when the ratio rx/vy = 1 and large beam size-rms mis- 
matched that the coupling produces an observable effect 
in the beam as a whole. This effect arises from both a 
beating in amplitude between the two coordinate direc- 
tions for the single-particle motion and from the core-core 
resonances, resulting in growth and transfer of emittance 
from one phase plane to another, in the beam developing 
an elliptical shape along directions preferential and there- 
fore in the fast transition from isotropic to anisotropic 
beam profiles. 

Nonlinear space-charge forces lead to equipartition- 
ing of energy between degrees of freedom. In space- 
charge dominated beams. Coulomb collisions are infre- 
quent to account for the energy transfer, whereas space 
charge waves have been shown to be a possible cou- 



pling canditate [l^, [13, H^. The equipartitioning of 
anisotropic beam involves nonlinear energy transfer and 
evolution towards a metaequilibrium state, as conse- 
quence of resonant phase mixing [4^ [g^, [Ml • We used 
FFTW (The Fastest Fourier Transform in the West) in 
order to compute frequencies numerically. The breath- 
ing mode frequency is given by ojx^ — 2.01334 and the 
quadrupole mode frequency is toxa — 2.01334. Both 
modes have the same frequency, therefore breathing 
mode and quadrupole mode are resonants. This core- 
core resonance are making the beam to develop an el- 
liptical shape with a increase in its size along the x- 
dircction [l5| . We computed numerically the frequencies 
ijjx and ujy of 2500 test particles. These 2500 test par- 
ticles have LOx = 1.05461 and 2455 test particles have 
ijjy — 1.05461. Thus 2455 test particles obey resonance 
condition Ux—t^y ~ 0. This number large of test particles 
in resonance causes emittance coupling. As shown in the 
histogram in Fig. 01 45 test particles have Uy = 1.00667 
. The particle-core (2:1) resonances are present in both 
quadrupole and breathing mode of the core oscillation, 
but only in y-direction of the test particles [2l|. Thus, 
halo is formed along y-direction. The anisotropy leading 
to coupling resonance [itI [isj in the presence of nonlin- 
ear space-charge forces was suggested as an approach to 
the equipartitioning question [2^ [13, [3l|, [s^ . 

The equipartition of the beam is driven for anisotropic 
processes. Based on turbulent equipartiton [s^ we pro- 
pose the anisotropic equipartition. Anisotropic equipar- 
tition corresponds to a phase space density uniform on 
the surface invariant of the (^ = {vyexY / {rx^yY) , where 
^ is the ratio of oscillations energies in the x and y di- 
rections, a version of the ergodic hypothesis where the ^ 
invariant play the role of the conserved energy [6^ . In 
the state of anisotropic equipartition, the temperature is 
stationary, the entropy grows in the cascade form, there 
is a coupling of transversal emittance, the beam develops 
an elliptical shape with a increase in its size along one 
direction and there is halo formation along one direction 
preferential. 

In order to understand the kinetic and dynamical be- 
havior of an anisotropic beams, in particular, to study 
mechanisms of equipartition connected with phase space, 
we compute the variables (volume, area of surface, area 
projected, temperature, heat flux and beam entropy) 
which characterize the anisotropic beam in phase space. 
We observed transverse phase planes coupling in a; — a; 
and y — y area projected defined by ax and Uy. This 
coupling is characterized by the energy exchange between 
the directions. Energy/acceptance exchange requires res- 
onant coupling. Energy /acceptance flow is due to the 
correlation jj between the velocity and the position of 
the particle in regions where the beam size contracts or 
expands. When the beam radius expands, the flow com- 
ponent has an outward direction and the thermal com- 
ponent decreases, when the beam radius contracts, the 
flow component is inward and the thermal component 
increases. The propagation of the beam is thus charac- 
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terized by the variation in the envelope beam which cor- 
relates with the generation of the heat flow and the vari- 
ation of the beam temperature. The elementary mech- 
anism responsible for the heat transfer from one degree 
of freedom to another is constituted by the effect of in- 
teractions between the individual particles and waves of 
collective fields [soj . 

The temperature Teg oscillates after steady with small 
fluctuations and entropy S grows by cascades as shown in 
Fig.O We observe to have different regimes depending on 
the value of time scales. There is first a phase of violent 
relaxation leading to a quasi-stationary state. This phase 
is followed by a thermalization leading to the Teq station- 
ary due to the thermal bath (field waves) i.e., the com- 
bined effect of imposed friction and diffusion [8l|, [H, [s^ , 
the diffusion process are arising from the fluctuations of 
the self-fields and the friction results from a polariza- 
tion process; not to collisions effects. The particles dy- 
namic described by these processes has a complex phase 
space structure. The effectively accessible phase space 
of the beam have a complicated geometrical structure. 
We showed that the dynamics of the volume /i and the 
surface a are coupled. They are not constant in func- 
tion of the axial coordinate of a beam s. But the surface 
area and volume become stationaries in function of the 
anisotropy variable x- 

This series of phenomena suggest an algorithm for fu- 
ture research. First, to map the building blocks of the 
beam core in the phase space that energy equipartition 
of the beam. The dynamic is described as a trajectory 
in the phase space. If the system has some integrals of 
motion, the trajectory is confined in the subspace which 
conserves the integrals. The system is defined as er- 
godic if the time average of any dynamical quantity is 
equal to the spatial average over the subspace. This sub- 
space is referred to as the ergodic region. At the equi- 
librium, energy is equally distributed to all degrees of 
freedom, which is called equipartition. In particular, this 
anisotropic beam has a quasi-stationary state that seems 
to suggest the existence of additional integrals, which 
conserves the anisotropy [H, . The anisotropy lead to 
the equipartition beam. 

Second, to describe self-consistently, the interactions 
between the individual particles and waves of collective 
fields via hamiltonian formalism [s^ . these interactions 
form the halo. The traditional description of the in- 
teracting particles and fields rests on the coupled set 
of Vlasov-Poisson equations. In the self-consistent ap- 
proach to the wave-particle interaction complements this 
usual treatment in that, to capture the physical mech- 
anism at work, the beam is partitioned in two popula- 
tions: core and tail. The idea behind this discrimination 
is simple: wave-particle interaction involves the resonant 
tail particles whose velocity is close to the phase velocity 
of the wave under consideration. These waves are just 
the collective macroscopic degrees of freedom, captur- 
ing the oscillations of other nonresonant, core particles, 
so that these core particles participate in the effective 



wave-particle dynamics only through the waves. 

Third, to analyze the halo-core interaction via ther- 
mal bath, ie, the Fokker-Planck equation. The Lenard- 
Balescu equation can be used to describe the evolution 
of the halo particle in a bath of field particles (core) at 
equilibrium. In that case, we have to consider that the 
distribution of the bath is given, that is f{v,t) = fo{v) 
where fo{v) is a stable stationary solution of the Vlasov 
equation (bath distribution). This procedure transforms 
the Lenard-Balescu equation into a Fokker-Planck equa- 
tion for the density probability P{v, s) of finding the halo 
particle with velocity v at time s. Chavanis [65] studied 
the relaxation of a test particle in a bath of field particles. 
Specifically, he considered a collection of N particles at 
statistical equilibrium (thermal bath) and introduced a 
new particle in the system. These results are well-known 
when the system is spatially homogeneous and memory 
effects can be neglected, as in the case of plasma physics. 
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